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Abstract 

Predictions for angular correlations between an arbitrary number of partons 
are derived in the high energy limit. The quantities considered depend on angles 
and primary energy through a single variable e which implies certain scaling 
properties and relations between quite different observables. These asymptotic 
predictions derived in the double log approximation of QCD are checked against 
Monte Carlo calculations at the parton and hadron level. 

INTRODUCTION 

Many new results from the detailed study of multiparticle production processes 
have been obtained in recent years, still there is no fully satisfactory model to describe 
all these phenomena. One approach is based on perturbative QCD, particularly suited 
for hard processes involving large momentum transfers, in which an initially scattered 
parton evolves by gluon Bremsstrahlung into a jet of partons and finally into the ob- 
servable hadrons, whereby the effect of the color confinement force is modeled. Another 
approach relies on a statistical and thermodynamic description in which the microscopic 
degrees of freedom are integrated over and only global quantities are kept. Hadrons 
are produced from a quark gluon plasma after a phase transition. This approach is 
applied in particular to nuclear collisions in the search for the quark gluon plasma. 
As many aspects of particle production change only moderately when going from the 
more elementary collisions like e + e~ annihilations to the more complex pp or nuclear 
collisions it is desirable to develop both the perturbative and the statistical methods 
and to explore their predictive power. 

In any case, for a satisfactory model of multiparticle production we would like to 
have finally a mathematical model based on a simple principle, there should be only 
a few parameters, some results should be obtained analytically - even if only approxi- 
mately - otherwise the structural properties of the theory can hardly be fully explored. 
A nice example of a model of this type is Hagedorn's bootstrap model: it is based 
on a single guiding principle, the bootstrap principle, the only parameters being the 



particle masses and the interaction volume; analytical results on multiparticle observ- 
ables can be derived from an equation for the generating functional of multiparticle 
densities.0 Although it has been attempted to extend the statistical approach to hard 
processes! the most convincing results for such processes are obtained today from the 
parton cascade model as derived from perturbative QCD. 

In this talk I would like to report on results from perturbative QCD on angu- 
lar correlations inside a parton jet which are obtained in collaboration with Jacek 
Wosiek.i'S The essential parameters of the parton cascade are the scale parameter A, 
which determines the coupling strength and a cutoff Q D which regulates the infrared 
and collinear singularities of the gluon Bremsstrahlung. As to the hadronization process 
we follow here the idea that it is sufficiently soft and the distribution of hadrons follows 
largely the distribution of partons.1 In particular for momentum spectra of partons and 
hadrons such an equality up to a constant has been established, if the parton cascade 
is evolved down to the hadronic (i.e. pion) mass scale Q Q pa Such a similarity may 
also be expected if the observable considered is "infrared save", i.e. does not depend 
explicitly on the cutoff Q Q . It is one of the motivations of the present study, as to what 
extent such a hypothesis of soft hadronization is actually correct and supported by the 
experiment. It would allow for a certain class of observables to be calculated just only 
in terms of Q Q and A parameters, resulting in a scheme of high predictive power. 



THEORETICAL SCHEME 

Our calculations are based on the double logarithmic approximation (DLA) of 
QCD. The probability to emit a single gluon of momentum K from a primary parton 
a of momentum P is given by 

ka (zr\jZif 2/7^. \dK dQpK d$ PK m 
Mp ta {K)d K = c a a {KQ PK ) — — — , (1) 
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where a 2 = r j 2 = Qa s /it is the QCD anomalous dimension of the multiplicity evolu- 
tion and c g = 1 or c q = 4/9 for gluon or quark jets respectively. For the running 
coupling constant we write a 2 (pr) = P 2 / '( m (pr /Qo) + A) whereby A = ln(Qo/A) and 
(3 2 = 12(^^c — §^/) _1 ~ 1.565 for 5 flavours Nf. In this approximation the recoil 
effects are neglected (i.e. energy and momentum conservation are violated), integrals 
are performed by retaining only the leading contributions from collinear and soft di- 
vergencies. The interference of the soft gluons is taken into account by the "angular 
ordering" prescription. In DLA one obtains the leading asymptotic behaviour, but the 
non-asyptotic corrections of relative order *JoT s are potentially large. 

The generating functional which includes the "leading" primary parton in the 
final state is given by the non-linear integral equationl 



Z P , a {u) = u{P) exp f [ M P , a (K)[u(K)Z K Ju} - l]d 3 K 



(2) 



In this form the virtual corrections are included which ensures the proper normalization 
Zp{u} = 1. The density distribution of n partons is then obtained by functional 

u=l 

differentiation after the test functions u(k) 

p^ih, k n ) = 5 n Z{u}/Su(k 1 )...Su(k n ) \u=l ■ (3) 



Similarly, the cumulant (connected) correlation function is derived as in Eq.(|3|) but with 
Z replaced by InZ. Starting with n = 1 one obtains from (^) and (fj) linear integral 



equations which can be solved recursively. Here we specialize on angular distributions 
and the corresponding, equations are obtained after integration over the momenta. We 
have obtained results!' @ on the general inclusive cumulant correlation function of n 
particles in their spherical angles. Of special interest are the distribution in the relative 
polar angle $12 between two partons and the multiplicity moments of general order n 
for particles falling into sidewise angular regions. 

It turns out that all these angular observables can be constructed from a generic 
function (5, P) which fulfils the integral equation 

hM(6,0,P) = dF>{6,0,P)+ f P ^ r^-a 2 (p T )h^(8,^,K). (4) 
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where 5 and d are the small and the large angles of the respective problem and the 



inhomogeneous term behaves at high energies like d^ ~ exp(2n/3y ln(P<5/A)). One 
finds that the natural variables of the problem, instead of P, d are 



e = ln(0/<$)/ ln(Ptf/A), C = l/(/3^/ln(Pi?/A)). (5) 

As ■d > 5 > Qo/P > A/P, we find < e < 1. The solution for In can be written as 
an expansion in £ ~ v/aj- In the high energy limit (e fixed, P — > 00) one obtains 

h in \5, P) ~ exp(2/3 v /ln(P^/A)cj(e, n)). (6) 

The scaling function u(e,n) is known in implicit form.! For small e it has a power 
expansion u(e, n) = n — (n 2 — l)e/2n + . . . Another useful approximation obtains from 
an expansion in 1/n which yields u(e,n) ~ ny/1 — e (1 — ln(l — e)) and has a 1% 
accuracy for e < 0.95 already for n = 2. An asymptotic behaviour of type (JD was also 
found in the study of azimuthal particle correlations.! 



POLAR ANGLE CORRELATIONS 

First we discuss the correlations p^($i2, P, ©) = dn/dtd 12 in the relative polar 
angle $12 of two partons both inside a forward cone around the initial parton of half 
opening angle O. For the study of scaling properties it is more convenient to consider 
the distribution in the variable e = ln(O/$i2)/ln(P0/A) given by p^{e) = dn/de = 
$12 ln(PG/A)/o' 2 '(i?i2). For the correlation f(e) = p^(e)/n 2 normalized by the multi- 
plicity in the coneB n one obtains in the high energy limit (e fixed, P — > 00) from Eq. 
dH) with h^ 2 '^^) ~ #12/? ($12) f° r either quark or gluon jet 

r(e) = 2/3 v /ln(PO/A)exp (-2/3 V / ln(PO/A) (2 - iv(e, 2))) . (7) 

Differences between quark and gluon jets and the influence of the leading particle show 
up at finite energies where we obtain 

f a (e) = c- 1 2 /exp(-y(2-^(£,2))-(c- 1 -l)yexp(-2y(l- >/T=I)) 

+ 2/3^/(^/(1 -e)J)exp(-y(2-vT^i)) (8) 

with ?/ = 2/3 v /ln(P0/A), / = 2pK (2pVX)/^r « 0.145 for A = In 2. To explore the 
scaling properties of r(e) we consider the quantity — lnf(e)/(2y / ln(PO/A)) depending 

*this normalization has better scaling properties! than the differential one used earlieJ] 



only on e for any momentum P or cone opening angle 0. The asymptotic limit and a 
finite energy result are shown in Fig. 1. As a test of our analytical calculations we also 
compared to the Monte Carlo evaluation of the parton cascade.! As can be seen from 
the figure, the predicted asymptotic scaling behaviour is nicely reproduced for small 
e < 0.5 in a large energy range (P > 20 GeV) whereas deviations occur for larger e 
(smaller relative angles i?i 2 approaching the cut-off Q /P). It should be noted that 
the normalization of the above expressions in (0, ||) and in particular the quark and 
gluon difference at finite energy are of non-leading order in the DLA and therefore less 
reliable. We have therefore adjusted the overall normalization of the functions in (|7|, 
H) to the Monte Carlo data. We also studied the effect of hadronization provided by 
the HERWIG MC. Again these effects are negligable for small e < 0.5.Q Preliminary 
data from DELPHI0 have recently given first evidence for approximate e-scaling in 
(0 > 30°) and a close proximity of the data to the parton model results. 

MULTIPLICITY MOMENTS FOR SIDEWISE ANGULAR 
CONE OR RING 

In a second application we consider particle multiplicities in a sidewise cone Sfl 
of half opening 5 at polar angle t? with respect to P and in an angular ring of width 
25 symmetrically around the primary parton direction, centered at polar angle *&. We 
will refer to the ring and the cone as the ID and 2D configurations. We calculate 
the factorial and cumulant multiplicity moments and c- n > or normalized by the 
multiplicities n in the respective angular regions F^ n ' = f^/n n and C^ 1 ' = c^/n n 
(j( 2 ) =< n{n — 1) >, = F^ — 1, etc.). The cumulant moments can be derived 
again from the generic equation (|J). For the normalized moments we find 

C {n) S) ~ ^/5) D{n - 1) exp (-2Py/ln(P&/A)(n - uj{e, n))) (9) 

with e = ln(??/5)/ ln(P??/A). The dependence on e is as in f(e) discussed above. Similar 
results have also been found by other groups.0lll In the modified LLA including next 
to leading effects in the exponent^ the power changes typically by 10%. 

The scaling properties of the various moments can again be conveniently investi- 
gated by projecting out the e-dependence of the exponent in Eq.@, one finds for 

M _ HW^"-^'] „( E , B) 

n 



n^ln(P^/A) 

in the high energy limit, or —C^ ~ 2/5(1 — y/1 — e) for large n independent of n,D. 

In Fig. 2a we plot — for the ring (D = 1) vs. s for different primary momenta 
P for the parton MC. There is a violation of e-scaling for small e but the asymptotic 
prediction is approached for high energies. We have repeated the same calculation as in 
(|ToD , but for the factorial moments F^ replacing in (|10D, see Fig. 2b. As F^ = 
+ 1 they approach the same asymptotic limit. Apparently, the nonasymptotic 
corrections are such that the scaling in e for small e sets in already at low energies for 
the factorial moments. The results for higher moments follow roughly the expectation 
([TOD , on the other hand the D = 2 moments show a more gradual dependence on e 
than predicted.^ 

^we used the program HERWIG§ with parameters A = 0.15 GeV, m q — m g = 0.32 GeV and 
without non-perturbative gluon splitting for the process e + e~ — > uu. 
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Figure 1: Rescaled polar angle correlation function for the high energy limit, Eq. (7) (full 
line), and for a quark jet of 45 GeV, Eq. (8) (dashed line), with normalization adjusted to 
the data as obtained from the HERWIG MC at the parton level. The insert shows the energy 
dependence of the same quantity for fixed e. 

An interesting aspect of these results is the remarkable universality of the various 
moments (different n, D) and also of the very different observable f(e) which all converge 
against the same limiting function after appropriate rescaling (see Figs. 1,2). 

SUMMARY 

The angular observables considered here, after appropriate rescaling, approach 
a limit in the normalized logarithmic angular variable e = ln(i?/5)/ln(Pi?/A) for e 
fixed, P — > oo. The comparison with the parton MC shows that this limit is already 
approached for 2-particle correlation r(e) and the factorial moments at present 
energies sufficiently far away from the angular cutoff Qo/P (e < 0.5) whereas the 
cumulant moments approach the asymptotic limit only at higher energies (P ~ 1 TeV). 
In the region of small e - not further discussed here - the observables approach a power 
behaviour in the angular variables and become independent of the cutoff Qo ("infrared 
safe"). It is in this region that also the hadronisation corrections are found small. 

The experimental confirmation of the universal high energy behaviour of rather 
different angular observables and of e-scaling with two redundant variables P and i? can 
provide further evidence for a soft confinement mechanism and parton hadron duality 
which allows to calculate hadronic distributions directly from the partonic ones. 
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Figure 2: (a) Rescaled cumulant moments for the ring as denned in Eq. (10) for the parton 
MC for different jet momenta P in comparison with the asymptotic prediction Eq. (10); 
(b) same as (a) but for factorial moments. 
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